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An extensive ab initio study of the physical properties of both linear and zigzag atomic chains 
of all Ad and 5d transition metals (TM) within the generalized gradient approximation by using 
the accurate projector-augmented wave method, has been carried out. The atomic structures of 
equilibrium and metastable states were theoretically determined. All the TM linear chains are 
found to be unstable against the corresponding zigzag structures. All the TM chains, except Nb, Ag 
and La, have a stable (or metastable) magnetic state in either the linear or zigzag or both structures. 
Magnetic states appear also in the sufficiently stretched Nb and La linear chains and in the largely 
compressed Y and La chains. The spin magnetic moments in the Mo, Tc, Ru, Rh, W, Re chains 
could be large (>1.0 /^a/atom). Structural transformation from the linear to zigzag chains could 
suppress the magnetism already in the linear chain, induce the magnetism in the zigzag structure, 
and also cause a change of the magnetic state (ferromagnetic to antiferroamgetic or vice verse). The 
calculations including the spin-orbit coupling reveal that the orbital moments in the Zr, Tc, Ru, 
Rh, Pd, Hf, Ta, W, Re, Os, Ir and Pt chains could be rather large (>0.1 /is/atom). Importantly, 
large magnetic anisotropy energy (>1.0 meV/atom) is found in most of the magnetic TM chains, 
suggesting that these nanowires could have fascinating applications in ultrahigh density magnetic 
memories and hard disks. In particular, giant magnetic anisotropy energy (>10.0 meV/atom) could 
appear in the Ru, Re, Rh, and Ir chains. Furthermore, the magnetic anisotropy energy in several 
elongated linear chains could be as large as 40.0 meV/atom. A spin-reorientation transition occurs in 
the Ru, Ir, Ta, Zr, La and Zr, Ru, La, Ta and Ir linear chains when they are elongated. Remarkably, 
all the 5d as well as Tc and Pd chains show the colossal magnetic anisotropy (i.e., it is impossible 
to rotate magnetization into certain directions). Finally, the electronic band structure and density 
of states of the nanowires have also been calculated in order to understand the electronic origin of 
the large magnetic anisotropy and orbital magnetic moment as well as to estimate the conduction 
electron spin polarization. 

PACS numbers: 73.63.Nm, 75.30.Gw, 75.75.+a, 61.46.-w 



INTRODUCTION 



Magnetism in nanostructures has been a very active 
research area in the last decades^— , because of its novel 
fundamental physics and fascinating potential applica- 
tions. Experimentally, modern methods of preparing 
metal nanowires have made it possible to investigate the 
influence of dimensionality on the magnetic properties. 
For example, Gambardella, et aiA : recently succeeded in 
preparing a high density of parallel atomic chains along 
steps by growing Co on a high-purity Pt (997) vicinal 
surface and also observed one-dimensional (ID) mag- 
netism in a narrow temperature range of 10~20 K. In 
the mean time, Li, et al£ reported that Fe stripes on the 
stepped Pd(110) substrate have a different magnetic easy 
axis than previous results. Structurally stable nanowires 
can also be grown inside tubular structures, such as the 
Ag nanowires of micrometer lengths grown inside self- 
assembled organic (calix[4]hydroquinone) nanotubes 6 . 
Short suspended nanowires have been produced by driv- 
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ing the tip of scanning tunneling microscope into contact 
with a metallic surface and subsequent retraction, leading 
to the extrusion of a limited number of atoms from either 
tip or substrate 7 -. Monostrand nanowires of Co and Pd 
have also been prepared in mechanical break junctions, 
and full spin-polarized conductance was observed 8 -. 

Theoretically, a great deal of research has been done 
on both finite and infinite chains of metal atoms. Theo- 
retical calculations at either semi-empirical tight-binding 
or ab initio density functional theory level for many in- 
finite/finite chains, e.g., linear chains of Co£~— , Fo n ' 14 , 
Ni, Pd±^, Pt, Cu&, AgiL±£, and Au^^-2 2 ., as well as 
zigzag chains of Fe^, Zr— and Av>2i, have been reported. 
Early studies of infinite linear chains of A u 21 ' 22 ' 24 ' 25 , 
Cu 18 , and Pd 20 have shown a wide variety of stable and 
metastable structures. Recently, the magnetic properties 
of transition metal infinite linear chains of Fe, Co, Ni, 
have been calculate d 10 ' 11 ' 13 ' 14 ' 26 . Possible magnetism in 
s- and sp-electron element linear and zigzag chains have 
also been studied theoretically. 27 These calculations show 
that the metallic and magnetic nanowires may become 
important for electronic/optoelectronic devices, quantum 
devices, magnetic storage, nanoprobes and spintronics. 

Despite of the above mentioned intensive theoretical 
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and experimental research, current understanding on the 
intriguing magnetic properties of nanowires and how 
magnetism depends their structural property is still in- 
complete. The purpose of the present work is to make a 
systematic ab initio study of the magnetic, electronic and 
structural properties of linear and zigzag atomic chains 
(Fig. 1) of all 4d and 5d transition metals (TM). Transi- 
tion metals, because of their partly filled d orbitals, have 
a strong tendency to magnetize. Nonetheless, only 3d 
transition metals (Cr, Mn, Fe, Co, and Ni) exhibit mag- 
netism in their bulk structures. It is, therefore, of interest 
to investigate possible ferromagnetic (FM) and antifcrro- 
magnetic (AF) magnetization in the linear chains of all 
Ad and 5d transition metals including Y, Zr, Nb, La, Hf 
and Ta zigzag chain which appear not to have been con- 
sidered. As mentioned before, recent ab initio calcula- 
tions indicate that the zigzag chain structure of, at least, 
Zr— , Ir— , Pti& and A u 16 ' 25 , is energetically more favor- 
able than the linear chain structure. Thus, we also study 
the structural, electronic and magnetic properties of all 
Ad and 5d transition metal zigzag chains in order to un- 
derstand how the physical properties of the monoatomic 
chains evolve as their structures change from the linear 
to zigzag chain. 

Relativistic electron spin-orbit coupling (SOC) is the 
fundamental cause of the orbital magnetization and 
also the magnetocrystalline anisotropy energy (MAE) of 
solids. The MAE of a magnetic solid is the difference 
in total electronic energy between two magnetization di- 
rections, or the energy required to rotate the magne- 
tization from one direction to another. It determines 
whether a magnet is a hard or soft one. Furthermore, 
it acts to reduce the magnitude of superparamagnetic 
fluctuation in nanostructures, and hence is a key factor 
that would determine whether the nanowires have po- 
tential applications in, e.g., high-density recording and 
magnetic memory devices. Ab initio calculations of the 
MAE have been performed for mainly the Fe and Co lin- 
ear chain a 10 ' 28 " — , while semiempirical tight-binding cal- 
culations have been reported for both linear chains and 
two-leg ladders of Fe and Co^Sr— . Very recently, we have 
carried out systematic ab initio calculations of both the 
MAE and also the magnetic dipolar (shape) anisotropy 
energy for all 3d transition metals in both the linear and 
zigzag structures^ Remarkably, although the SOC is 
rather weak in 3d transition metals, compared with Ad 
and 5d transition metals, we found that the FM Ni linear 
chain has a gigantic MAE of ^12 meV/atom^ There- 
fore, as a continuing endeavor to find nanowires with a 
large MAE, we have calculated the MAE and the mag- 
netic dipolar (shape) anisotropy energy for all Ad and 5d 
transition metals in both the linear and zigzag structures. 
Although in this paper we study only free-standing Ad 
and 5d transition metal chains, the underlying physical 
trends found may also hold for monoatomic nanowires 
created transiently in break junctions^ or encapsulated 
inside ID nanotubes^i 2 ^ or deposited on weakly inter- 
acting substrates 3 ^, albeit, with the actual values of the 



physical quantities being modified. 

The rest of this paper is organized as follows. In the 
next section, we briefly describe the theory and compu- 
tational details we used. The calculated structural and 
magnetic properties as well as band structures of the lin- 
ear Ad and 5d transition metal chains are presented in 
Sec. III. The calculated structural, magnetic and elec- 
tronic properties of the zigzag Ad and 5d transition metal 
chains in both equilibrium and local energy minimum 
states are reported in Sec. IV. The relative stability of 
the linear and zigzag chain structures is analyzed in Sec. 
V. The calculated magnetic anisotropy energies and mo- 
ments of both linear and zigzag chains are presented, 
and also discussed in terms of the calculated d-orbital- 
decomposed DOSs in Sec. VI. Finally, a summary is 
given in Sec. VII. 



II. THEORY AND COMPUTATIONAL 
METHOD 

In the present calculations, we use the accurate frozen- 
core full-potential projector augmented-wave (PAW) 
method^ as implemented in the Vienna ab initio simu- 
lation package (VASP)— The calculations are based 
on density functional theory with the generalized gradi- 
ent approximation (GGA)2&. The free-standing atomic 
chains are modelled by a two-dimensional array of infi- 
nite long, straight or zigzag wires. For both linear and 
zigzag chains, the nearest wire- wire distance between the 
neighboring chains is, at least, 15 A, which should be 
wide enough to decouple the neighboring wires. A large 
plane-wave cutoff energy of ~350 eV is used for all Ad 
and 5d transition metal chains. 

The equilibrium bond length (lattice constant) of the 
linear atomic chains in the nonmagnetic (NM), ferromag- 
netic (FM) and antiferromagnetic (AF) states is deter- 
mined by locating the minimum in the calculated total 
energy as a function of the interatomic distance. The 
results are also compared with that obtained by struc- 
tural optimizations, and the differences are small (within 
0.4 %) for, e.g., the Ru, Rh and Pd chains. For the 
zigzag chains, the theoretical atomic structure is deter- 
mined by structural relaxations using the conjugate gra- 
dient method. The equilibrium structure is obtained 
when all the forces acting on the atoms and the axial 
stress are less than 0.02 eV/A and 2.0 kBar, respectively. 
The T-centered Monkhorst-Pack scheme with a /c-mesh 
of 1 X 1 X n (n = 40) in the full Brillouin zone (BZ), in 
conjunction with the Fermi-Dirac-smearing method with 
cr = 0.01 eV, is used to generate Appoints for the BZ in- 
tegration. With this fc-point mesh, the total energy is 
found to converge to within 10~ 3 eV. 

Because of its smallness, ab initio calculation of the 
MAE is computationally very demanding and needs to 
be carefully carried out (see, e.g., Refs. |39jj40j). A very 
fine fc-point mesh with n being 200 for both the linear and 
zigzag chains, is used. The same fc-point mesh is used for 



FIG. 1: (color online) Schematic structure diagram for (a) the 
linear and (b) zigzag atomic chains. 



the band structure and density of states calculations. As 
in our previous publication on the 3d TM chains^, we 
use the force theorem approach to calculate the MAE, 
i.e., the MAE is calculated as the total energy difference 
between the two relativistic band structure calculations 
for the two different magnetization directions (e.g., par- 
allel and perpendicular to the chain) concerned using the 
frozen charge density obtained in a prior self-consistent 
scalar relativistic calculation. 41 The total energy conver- 
gence criteria is 10 -7 eV/atom. 



III. LINEAR CHAINS 
A. Magnetic state and spin magnetic moment 

The calculated equilibrium bond lengths (d) and spin 
magnetic moments of all the 4d and bd transition metal 
linear chains in the NM, FM and AF states are displayed 
in Fig. 2, and Fig. 3, respectively. They are also listed 
in Table I. The calculated total energy relative to that of 
the NM state (i.e., the magnetization energy) of the FM 
and AF linear atomic chains are also shown in Fig. 2, 
Fig. 3 and Table I. It is clear from Figs. 2 and 3 that all 
of the Ad and 5d TM elements except Y, Nb, La, Ta, Os 
and Pt, become magnetic in the linear chain structure. 
Furthermore, for all the 4d and 5d TM elements, except 
Y, Nb, La, Ta, Os and Pt, NM state is unstable and the 
ground state is either FM and AF (see Fig. 2, Fig. 3 and 
Table I). Among the 4d TM linear chains, the ground 
state for the Zr, Ru, Rh, and Pd chains is ferromagnetic 
while that for the Mo, and Tc chains is antiferromag- 
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FIG. 2: (color online) (a) Equilibrium bond length (A), (b) 
magnetization energy (AE) (i.e., the total energy of a mag- 
netic state relative to that of nonmagnetic state) (AE — 
e fm(af) _ E NM} and ( c ) spm magnetic moments (jits) of 
all the Ad TM linear atomic chains in the NM, FM, and AF 
states. 
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FIG. 3: (color online) (a) Equilibrium bond length (A), (b) 
magnetization energy (AE) (i.e., the total energy of a mag- 
netic state relative to that of nonmagnetic state) (AE — 
e fm(af) _ E NM} and ( c ) spm magnetic moments (jits) of 
all the 5d TM linear atomic chains in the NM, FM, and AF 
states. 



4 



netic. For the 5c? TM linear chains, the ground state 
for the Hf and Ir chains is ferromagnetic and the ground 
state for the Re and W are antiferromagnetic. The Y, 
Nb, La, Ta, Os and Pt chains are nonmagnetic at the 
equilibrium bond length. We recently reporteo-2, that in 
the 3d TM linear chains, the equilibrium bond length in 
a magnetic state is significantly larger than that in the 
nonmagnetic state. For example, the magnetization in- 
duced increase in the bond length in the Cr chain is 54 
%. In contrast, Fig. 2, Fig. 3 and Table I show that in 
the Ad and 5d TM linear chains, the difference in bond 
length between a magnetic (FM or AF) state and the NM 
state is much smaller. The largest lattice magnetolattice 
expansion occurs in the AF Re chain but it amounts only 
to 3 %. This is due to much weak magnetization in the 
Ad and 5d linear chains, as indicated by the smaller mag- 
netic moments and much smaller magnetization energies 
in these atomic chains at equilibrium (Table I). 

To see how the magnetic properties of the atomic 
chains evolve with the interatomic distance, we plot the 
spin and orbital moments for some Ad (Y, Zr, Nb, and 
Pd) and 5d (La, Hf, Ta, Os, Ir and Pt) TM chains in 
the FM state as a function of the bond length in Fig. 
4. For most selected ferromagnetic TM (except Y, La, 
Zr and Hf) chains, the spin moment generally becomes 
larger as the bond length is increased from its equilib- 
rium value (Table I). Interestingly, the spin moment of 
the Hf chain, in contrast, decreases monotonically when 
the chain is elongated, and eventually disappears at the 
bond length of 2.7 A [Fig. 4(c)]. The spin moment of the 
Zr chain decreases slightly as the bond length increases, 
but increases again when the bond length goes beyond 
~2.7 A. Surprisingly, when the Y (La) chain is sufficiently 
compressed [at the bond length of ^2.35 (2.65) A], the 
ferromagnetism appears, and the spin moment increases 
as the chain is further compressed. Finally, for the Y, 
La, Ta, Os and Pt chains, the sufficient elongation of the 
bond length would induce a FM state (Fig. 4). 

Our calculated bond lengths, spin moments and mag- 
netization energies generally agree rather well with 
available previous ab initio calculation a 15 ' 20 i 26 ' 29 i 42 ~ — . 
Nonetheless, a few notable differences exist. For example, 
our calculated bond length (2.13 A) of the AF Mo chain is 
7.0 % smaller than that (2.28 A) reported in Ref. but 
in good agreement with Ref. [H (2.15 A). Also, our equi- 
librium bond lengths of the Os and Pt chains are smaller 
than that from Ref 20;, but the differences are within 2.7 
%. Another notable difference is that our calculations 
suggest that the Os chain is nonmagnetic in equilibrium 
but become ferromagnetic only when the bond length is 
larger than ^2.55A (Fig. 4), while, according to Ref. [2p| . 
it is ferromagnetic at the equilibrium bond length. Our 
calculated magnetization energies (Table I) for the AF 
Mo and Tc chains are smaller than that reported in Ref. 
(197 and 53 meV/atom, respectively) and in Ref. [l5| 
(92 and 65 meV/atom, respectively). 



TABLE I: Equilibrium bond lengths (d) (in A), total energies 
(Et) (in meV/atom) in the FM and AF states (relative to the 
NM state), and spin magnetic moments (m s ) (in /is/atom), 
of the Ad and 5d transition metal linear chains from scalar 
rclativistic calculations. 





dNM 






dFM 






dAF 










Ad metals 
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2.95 














Zr 


2.54 


-3.29 


0.628 


2.54 








Nb 


2.34 














Mo 


2.09 








-69.53 


1.317 


2.13 


Tc 


2.19 








-21.02 


1.268 


2.23 


Ru 


2.21 


-28.30 


1.118 


2.25 








Rh 


2.25 


-9.19 


0.328 


2.25 








Pd 


2.43 


-0.05 


0.684 


2.46 








Ae 


2.66 






5d metals 








La 


2.98 














Hf 


2.60 


-0.46 


0.137 


2.60 








Ta 


2.40 














W 


2.29 








-22.59 


1.465 


2.34 


Re 


2.26 








-118.52 


1.729 


2.32 


Os 


2.25 














Ir 


2.28 


-27.32 


0.660 


2.28 








Pt 


2.38 














Au 


2.60 















TABLE II: Spin (m s ) and orbital (m ) magnetic moments (in 
fiB /atom) of the magnetic Ad and 5d transition metal linear 
chains at the equilibrium bond lengths (Table I) with mag- 
netization parallel (m|| z) and perpendicular (ml z) to the 
chain axis from fully relativistic charge selfconsistent calcula- 
tions. 

mil z m_L I 







m s 


m D 


m s 


TO 






Ad 


metals 






Zr 


(FM) 


0.631 


-0.065 


0.610 


-0.007 


Mo 


(AF) 


1.337 


-0.008 


1.181 


0.005 


Tc 


(AF) 


1.353 


0.463 


1.252 


0.046 


Ru 


(FM) 


1.115 


-0.106 


1.076 


0.058 


Rh 


(FM) 


0.317 


0.428 


0.017 


0.000 


Pd 


(FM) 


0.345 


-0.043 


0.636 


0.126 






5d 


metals 






Hf 


(FM) 


0.235 


-0.198 






W 


(AF) 


1.184 


-0.307 


1.371 


-0.005 


Re 


(AF) 


1.564 


0.115 


1.644 


0.146 


Os 


(FM) 






0.444 


0.046 


Pt 


2.38 


0.124 


0.100 







B. Orbital magnetic moment and colossal magnetic 
anisotropy 

The spin and orbital magnetic moments in the mag- 
netic Ad and 5d TM atomic chains in equilibrium from 
the fully relativistic charge selfconsistent calculations are 
listed in Table II. We note that the SOC affect slightly 
the spin moments in the AF TM chains and also FM Zr 
and Ru chains (see Tables I and II). However, in the other 



5 




2 2.2 2.4 2.6 2.8 
Interatomic distance (A) 



2.2 2.4 2.6 2.8 3 
Interatomic distance (A) 



FIG. 4: (color online) Spin (left panels) and orbital (right pan- 
els) magnetic moments as a function of interatomic distance 
of the ferromagnetic Y, La, Zr, Hf, Nb, Ta, Os, Ir Pd and Pt 
linear chains. In the left panels, "no-SOC" means the results 
from the scalar relativistic calculations. "Para" ("perp") de- 
notes the magnetization being parallel (perpendicular) to the 
chain axis. The spin magnetic moment for the Pd chain goes 
to zero at ~ 3.6 A. 



cases, the spin magnetic moments in Table II are gener- 
ally much smaller than that obtained from the scalar- 
relativistic calculations (Table I), unlike in the 3d TM 
chains where the SOC hardly affects the spin magnetic 
moments^, In fact, the SOC completely suppresses the 
spin magnetic moment in the FM Ir chain in equilib- 
rium (Tables I and II). Interestingly, the SOC- induced 
reduction of the spin magnetic moment is magnetization- 
direction dependent. Table II shows that the spin mo- 
ment of the Rh chain with magnetization parallel to the 
chain axis remains almost unchanged while that perpen- 
dicular to the axis becomes nearly diminished. In the 
Pd chain, in contrast, the spin moment for magnetiza- 
tion along the axis decreases nearly by half while that 
perpendicular to the axis remains nearly unchanged (Ta- 
ble II). Dramatically, in the FM Hf chain, the SOC fully 
suppresses the magnetization when the magnetization is 
perpendicular to the chain axis, but it nearly doubles the 



spin moment when the magnetization is along the axis. 
This interesting magnetic anisotropy is called the colos- 
sal magnetic anisotropy (CMA) by Smogunov et alM, 
who reported recently this CMA in the Pt chain. The 
CMA means that a magnetization magnitude could be 
finite only along certain directions and also that it is 
strictly impossible to rotate magnetization into certain 
directions. Earlier calculations^ also suggested the CMA 
to occur in the Rh linear chain. Our calculations here not 
only corroborate this finding of Smogunov et but 
also reveal the CMA in other 5c? transition metal linear 
chains such as Hf and Os (Table II). 

When the SOC is not taken into account, the spin mo- 
ment in the 4c? and 5d TM linear chains generally in- 
creases monotonically as the bond length is increased, as 
shown in Fig. 4. However, the behavior of the magnetic 
properties of the 5g? TM linear chains under the influ- 
ence of the SOC is very different from that of the 3d TM 
chains. For example, from the scalar relativistic calcu- 
lations, the Ir chain at the interatomic distance starting 
from 2.0 to 3.0 A, has a finite magnetic moment in the 
range of 0.4—2.4 /is/atom (Fig. 4g). When the SOC 
is taken into account, the Ir chain becomes nonmagnetic 
when the interatomic distance is smaller than 2.5 A, but 
has a finite magnetic moment when the interatomic dis- 
tance larger then 2.5 A. Similar behavior can also be 
seen in the Os, and Pt chains (Fig. 4 g and i). Our 
scalar relativistic calculations show that the Os chain is 
nonmagnetic if the interatomic distance is below 2.4 A, 
whilst our fully relativistic calculations indicate that, for 
the magnetization perpendicular to the chain direction, 
it become magnetic at the interatomic distance above 2.2 
A. In contrast, for the axial magnetization, the Os chain 
would become ferromagnetic only when the interatomic 
distance is larger than 2.5 A. Therefore, the Os chain 
exhibits the CMA^ when the interatomic distance falls 
between 2.2 and 2.5 A. Fig. 4 further shows that the 
Hf, Ir and Pt linear chains also exhibit the CMA in the 
interatomic distance of 2.25-3.0 A(Hf), 2.52-2.78A(Ir) 
and 2.30-2.63 A(Pt), respectively. 

The SOC provides the essential symmetry breaking 
that gives rise to orbital magnetization in magnetic 
solids. When the SOC is included in our calculations, 
the calculated orbital magnetic moments in the FM Ad 
TM chains at equilibrium bond length are listed in Table 
II. Surprisingly, even for the 5c? TM linear atomic chains, 
the calculated orbital magnetic moments are not large. 
For example, the calculated orbital moments in the 5d 
TM linear atomic chains are within — 0.2 /is/atom (Ta- 
ble II). The calculated orbital moments in the Ad TM 
linear chains in equilibrium can be larger, e.g., being — 
0.2 /x^/atom in the AF Tc and FM Rh chains with the 
magnetization along the chain direction (Table II). re- 
spectively Therefore, although the SOC is stronger in 4c? 
and 5c? transition metals than in 3c? ones, the calculated 
orbital magnetic moments in the 4c? and 5c? transition 
metals chains at the equilibrium bondlength is not neces- 
sarily larger than in the 3c? transition metal chains^ As 
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FIG. 5: Band structures of the Zr (upper panels) and Ir 
(lower panels) linear chains at 2.6 A. Left panels: the scalar- 
relativistic band structures; the middle and right panels: the 
fully relativistic band structures with the the magnetization 
parallel to and perpendicular to the chain axis, respectively. 
In the left panels, the solid and dashed lines represent (spin 
up) and (spin down) bands, respectively. The Fermi level (the 
dotted horizontal line) is at the zero energy. 

for the spin moments, the magnitude of the orbital mo- 
ments generally increases monotonically with the bond 
length, as can be seen in Fig. 4, with one notable ex- 
ception of the La chain (Fig. 4b). The orbital moment 
shows a strong dependence on the magnetization orien- 
tation (Fig. 4, right panels). As in the 3d TM chains^, 
the orbital moment in the 4d and 5d TM chains with the 
magnetization along the chain direction is usually much 
higher than that for the magnetization perpendicular to 
the chain. However, in the Pd chain, the orbital mo- 
ment with the magnetization along the chain direction is 
significantly larger than that for the magnetization per- 
pendicular to the chain. (Fig. 4j). 



C. Band structures and density of states 

Let us now examine the band structure of selected 
transition metal linear chains in order to understand the 
calculated magnetic properties. The energy bands ob- 
tained without and also with the SOC for the Zr and Ir 
chains in the FM state at 2.6 A are plotted in Fig. 5. In 
the absence of the SOC, because of the uniaxial rotational 
symmetry, the bands may be grouped into three sets, 



namely, the nondegenerate s- and d z 2 -dominant bands, 
double degenerate (d xz , d yz ), and (d x 2_ y 2, d xy ) dominant 
bands (see the left panels in Fig. 5). The {d x 2_ y 2, d xy ) 
bands are narrow because the d x 2_ y 2 and d xy orbitals are 
perpendicular to the chain, thus forming weak S bonds. 
The (d xz , dy Z ) bands, on the other hands, are more dis- 
persive due to the stronger overlap of the d xz and d yz 
orbitals along the chain, which gives rise to the tt bonds. 
The s- and d z 2 dominant bands are most dispersive since 
these orbitals form strong a bonds along the chain. The 
left panels in Fig. 5 show that the less dispersive d x 2_ y 2 ; 
d xy bands are near the Fermi level and spin-split. In the 
Zr linear chain, one spin-split d x 2_ y 2, d xy band is par- 
tially occupied near the T-point while the other band is 
completely empty. In the Ir chain, one split band lies 
completely below the Fermi level while the other band 
is only partially occupied. Thus, the relatively narrow 
d x 2_ y 2, d xy bands play an important role in magnetism, 
and that is the main reason why Zr and Ir chains are 
ferromagnetic at the bondlength of 2.6 A. 

The directional dependence of the orbital magnetiza- 
tion can be explained by analyzing the fully relativistic 
band structures (see Fig. 5). For the Zr linear chain 
with the axial magnetization (Fig. 5b), the doubly de- 
generate d x 2_ y 2, d xy bands are split into two with angu- 
lar momenta to/ = ±2. If one of them is fully occupied 
and the other is empty, the resulting orbital moment is 
2. Nonetheless, in the Zr linear chain, both are par- 
tially occupied with different occupation numbers (Fig. 
5b), resulting in an orbital moment of -0.07 /is/atom. Of 
course, the larger the SOC splitting, the larger the differ- 
ence in the occupation number and hence the larger the 
orbital moment. Therefore, the Ir chain has a larger axial 
orbital moment (1.06 /is /atom), because one of the split 
d x 2_ y 2 ; d xy bands lies almost completely below the Fermi 
level (see Fig. 5e). However, for the perpendicular mag- 
netization, the d x 2_ y 2, d xy bands remain degenerate (Fig. 
5c and Fig. 5f ) and hence do not contribute to the orbital 
magnetization. Nonetheless, as pointed out in Ref. H . 
the SO-split d x 2_ y 2, d xy -(d XZl d yz ) bands near the Fermi 
level would hybridize (Fig. 5f) and this hybridization 
would give rise to a smaller perpendicular orbital mo- 
ment of 0.27 /us/atom in the Ir linear chain. For the Zr 
chain, this hybridization does not occur near the Fermi 
energy (Fig. 5c). Therefore, the Zr chain have a tiny or- 
bital moment of -0.01 /is/atom when the magnetization 
is perpendicular to the chain axis. Of course, when the 
SOC is included, the degenerate d xz , d yz bands are also 
split into the to/ = — 1 and +1 bands for the axial mag- 
netization, but remain degenerate for the perpendicular 
magnetization (see Fig. 5). This SOC splitting of the 
{d xz , d yz ) band and {d x 2_ y 2 1 d xy ) band is proportional 
to | < d xz \H S o\d yz > | 2 and | < d x 2_ y 2 \H S o\d xy > | 2 , 
respectively. Here Hso is the SOC Hamiltonian. Since 
I < d xz \H so \d yz > | 2 :| < d x 2_ y 2\H SO \d xy > | 2 = 1:4, 46 
the SOC splitting of the (d xz , d yz ) bands is much smaller 
than the (d x 2_ y 2, d xy ) bands (see Fig. 5). Therefore, the 
(d xz , d yz ) bands would make a much smaller contribution 
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TABLE III: Numbers (nt and ni) of the spin-up and spin- 
down conduction bands crossing the Fermi level, and spin- 
polarization P at the Fermi level for the Ad and 5d TM atomic 
chains in the FM state. 





linear chain 


zigzag 


chain 




(nt, ni) 


P 


{nt, ni) 


P 






Ad metals 






Y 






(2,3) 


0.09 


Zr 


(4,3) 


0.48 


(3,4) 


-0.09 


Mo 








-0.11 


Ru 


(3,4) 


-0.41 






Rh 


(3,4) 


-0.65 


(3,6) 


-0.53 


Pd 


(3,4) 


-0.22 


(3,6) 


-0.30 






5d metals 






Hf 


(4,3) 


-0.31 






W 






(6,8) 


-0.16 


Re 






(2,4) 


-0.38 


Os 






(5,7) 


-0.10 


Ir 


(3,4) 


-0.47 


(7,9) 


-0.28 




-4 -2 
Energy (eV) 



4-202 
Energy (eV) 



to the orbital magnetization 

Electric and spin current transports are determined by 
the characteristics of the band structure near the Fermi 
level (Ep) in the systems concerned. Therefore, it would 
be interesting to examine the energy bands and density 
of states (DOS) of the atomic chains in the vicinity of 
the Ep. The spin-decomposed DOS for all the Ad and 
5d linear chains in equilibrium are displayed in Fig. 7 
and Fig. 8, respectively. For the FM Zr, Ru, Rh, Pd, Hf 
and Ir chains, the density of states at the Ep are spin- 
polarized (Fig. 7 and Fig. 8). This is usually quantified 
by the spin-polarization P defined as 



P = 



N t (E F ) - N^Ep) 
N t (EF) + N i {E F )'- 



(1) 



where N^(Ep) and N±(Ep) are the spin-up and spin- 
down DOS at the Ep, respectively. The most use- 
ful materials for the spintronic applications are the so- 
called half-metallic materials in which one spin channel 
is metallic and the other spin channel is insulating. The 
spin-polarization for these half-metals is either 1.0 or -1.0, 
and the electric conduction would be fully spin-polarized. 
The calculated spin-polarization and also the numbers of 
the conduction bands that cross the Fermi level in the 
Ad and 5d TM chains are listed in Table III. It is clear 
that the P of the FM Zr, Ru, Rh, and Ir linear chains 
is rather large (> 0.4), though still smaller than many 
3d TM linear chains^ None of the Ad and 5d TM lin- 
ear chains in the FM state is half-metallic. Interestingly, 
the FM Zr chain has a positive spin polarization, while 
the Ru, Rh, Hf, Pd and Ir chains have a negative spin 
polarization (Table II, Fig. 7 and Fig. 8). 



FIG. 6: (color online) Density of states of the FM Ad TM 
linear atomic chains at the equilibrium bond length. The 
Fermi level (dotted vertical lines) is at the zero energy. 




1 1 1 1 1 1 1 1 1 1 , 1 1 1 , 1 




- Ir I!) . 
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-4 -2 
Energy (eV) 



-4-2 2 
Energy (eV) 



FIG. 7: (color online) Density of states of the FM 5d TM 
linear atomic chains at the equilibrium bond length. The 
Fermi level (dotted vertical lines) is at the zero energy. 



IV. ZIGZAG CHAINS 

The zigzag structure for metal monoatomic wires has 
already been observed in experiments^. Among Ad 
and 5d transition metals, structura l 1 ^ 23 and magnetic^ 
properties of Zr, Rh, Pd, W, Ir, and Pt zigzag atomic 
chains have been studied theoretically in recent years. 
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TABLE IV: Equilibrium structural parameters (see Fig. lb 
for symbols di, d,2, a), spin magnetic moment (m s ) and mag- 
netization energy (AE) of the 4d and 5d transition metal 
zigzag chains from the scalar relativistic calculations, di and 
d2 are in the unit of A, and a is in the unit of degree. AE is 
in the unit of meV/atom, and m s in the unit of /^/atom. 
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388.1 




fNivn 

l IN 1V1 1 


9 81 


9 71 


7 

00. / 








1 r 1V1 1 


9 8^ 


9 7/1 


^8 fi 
Oo.O 


n 1 fi9 

U. IOZ 


9 


Nb 


(NM) 


2.51 


2.63 


61.5 






IVIO 




9 


z.Ol 


01. ( 
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2.40 


9 47 
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(FM) 


2.49 


2.48 
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9 4fi 


fift 7 
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Rh 


(NM) 


2.39 


2.60 


62.6 








fFivn 

1 r 1V1 1 
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9 /IQ 


oy. 1 


1.000 


-ou.u 


Pd 


(NM) 


2.56 


2.64 


61.0 








1 r ivi ) 




9 fifi 


fii ^ 


^Q9 


n 7 




(At ) 


9 ^fi 
Z.OD 


9 fi/1 


ah q 
ou.y 


n 9fifi 
u.zoo 


9 9 


As? 


(NM) 


2.73 


2.78 

5c? mcteils 


60.9 






La 


(NM) 


3.24 


0. 1U 


00.0 






Hf 




2.89 


2.71 


57.7 






Ta 


(NM) 


2.75 


2.51 


56.7 






W 


(NM) 


2.48 


2.56 


61.0 








(FM) 


2.48 


2.56 


61.0 


0.262 


-1.8 


Re 


(NM) 


3.22 


2.25 


48.8 








(FM) 


3.22 


2.25 


44.3 


0.516 


-259.5 


Os 


(NM) 


2.49 


2.44 


59.3 








(FM) 


2.50 


2.44 


59.1 


0.457 


-26.7 




(AF) 


2.50 


2.44 


59.2 


0.360 


-7.9 


Ir 


(NM) 


2.45 


2.53 


61.0 








(FM) 


2.44 


2.56 


61.5 


0.629 


-36.7 


Pt 


(NM) 


2.49 


2.65 


62.0 






Au 


(NM) 


2.67 


2.76 


61.1 







In the present paper, we perform a systematic ab initio 
study of the structural, electronic and magnetic proper- 
ties of the zigzag chain structure of all the Ad and 5d 
transition metals. 



A. Structure and magnetic moments 

The calculated equilibrium structural parameters (Fig. 
lb), spin magnetic moment and magnetization energy of 
the Ad and 5d TM zigzag chains are listed in Table IV. 
First of all, Table IV shows that all the zigzag chains 
except the Re one, look like planar equilateral triangle 
ribbens, i.e., the two bond lengths d\ and d 2 are similar 
and the angle a is close to 60° (Fig. lb). The equilib- 
rium bond lengths d\ and d 2 are generally a few percents 
larger than the bond length d of the corresponding linear 
chains (Table I). This is because the zigzag chains which 



TABLE V: Structural parameters (di, ffe, a) (see Fig. lb) of 
the zigzag chain at the second local energy minimum state, di 
and d2 are in the unit of A, and a is in the unit of degree. AE 
(meV/atom) is the energy difference between the second en- 
ergy minimum and the corresponding energy minimum listed 
in Table IV. The second local minimum state of the Zr and 
Ir chains only is ferromagnetic with a spin moment of 0.295 
and 0.285 /is /atom, respectively. The elements whose zigzag 
chains do not have the second energy minimum are not listed 
here 







di 


d 2 


a 


AE 








Ad metals 






Zr 


(FM) 


4.25 


2.44 


29.6 


1015 


Nb 


(NM) 


2.80 


2.45 


55.1 


-41.6 


Mo 


(NM) 


3.10 


2.24 


46.3 


-295.7 


Tc 


(NM) 


3.32 


2.20 


41.1 


-24.0 








5d metals 






W 


(NM) 


3.00 


2.31 


49.6 


-127.8 


Re 


(NM) 


3.20 


2.21 


43.7 


-102.6 


Ir 


(FM) 


4.00 


2.27 


28.4 


207.4 


Pt 


(NM) 


4.27 


2.37 


25.8 


412.0 


Au 


(NM) 


4.60 


2.55 


28.2 


330.6 



form planar equilateral triangle ribbens, have a higher 
coordination number (four) than that (two) of the linear 
chains. Similarly, all these bond lengths are shorter than 
their counterparts in the bulk structures. For example, 
the bond lengths for bcc Nb, bcc Mo, foe Rh, fee Pd, bec 
W, foe Ir and foe Pt are, 2.86, 2.73, 2.68, 2.75, 2.86, 2.72 
and 2.77 A, respectively.— 

Our calculated equilibrium structural parameters 
(di,d 2 ,a) agree reasonably well with available previous 
calculation o 16 ' 17 ' 23 . For example, Lin et al. reported 
d x = 2.86 A, d 2 = 2.74 A, a = 58.5° for the Zr zigzag 
chain, being consistent with our values in Table IV. Re- 
ported parameters d\, d 2 , and a for the W (2.44, 2.59, 
61.9), Os (2.48, 2.56, 61.1), Pt (2.58, 2.73, 61.9) and Au 
(2.64, 2.73, 61.13) chains (estimated from Figs. 3 and 5 
in Ref. Il6l ) are in rather good agreement with our re- 
sults in Table IV. One exception is the Ir chaini£ where 
di = 2.50 A, d 2 = 4.53 A, a = 74.0° differs substantially 
from the present results. Secondly, all the Ad and 5d 
TM zigzag chains except that of Nb, Tc, La, Hf, Ta and 
Pt, have magnetic solutions in the equilibrium structures 
(Table IV). Further, the Zr, Mo, Ru, Rh, W, Re, Os, and 
Ir zigzag chains are most stable in the FM state, whilst 
the ground state of the Pd zigzag chain is antiferromag- 
netic. For comparison, the ground state of the linear Mo, 
Tc, W and Re chain is antiferromagnetic (Table I). The 
FM Ru and Rh zigzag chains have a rather large spin 
moment of ~1.5 /is/atom, though the other magnetic 
zigzag chains generally have a small spin moment (< 1.0 
/tie/atom) (Table IV). Interestingly, the ground state of 
the Y zigzag chain is nonmagnetic, though it has a FM 
solution with a spin magnetic moment of ~0.5 /^s/atom. 
Note that none of the Ad and 5c? TMs is magnetic in their 
bulk structures in nature. Thirdly, for some Ad and 5c? 
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transition metals, the ground state magnetic configura- 
tion changes when the structure changes from the linear 
to zigzag chain. For example, the Tc and Hf elements are 
nonmagnetic in their equilibrium zigzag chain structures, 
though they are, respectively, antiferromagnetic and fer- 
romagnetic in their equilibrium linear chain structures 
(Table I). This is due to the increase in the coordina- 
tion number in the zigzag chains because most of them 
form a planar equilateral triangle ribben. On the other 
hand, the Y and Os elements become ferromagnetic in 
the zigzag chains even though they are nonmagnetic in 
the linear chains. Finally, the ground state of the Mo, W 
and Re chains changes from the AF state in the linear 
chain to the FM state in the zigzag structure. 

The Zr^, Ruil, Os, Au±£ Ir and Pt 16 ' 55 zigzag chains 
were reported to have a metastable non-triangular elon- 
gated zigzag structure with a being ~ 30.0°. The ex- 
istence of this second energy minimum elongated zigzag 
structure (a being ~ 30.0°) is believed to be crucial to 
the formation of a transient atomic chain in the break- 
junction experiments . 16 ! 55 To systematically study these 
possible elongated zigzag structures, we therefore further 
calculated the total energy as a function of the fixed lat- 
tice constant d± with d\ varying from 2.0 A to 6.0 A 
for all the Ad and 5d zigzag chains. The structural pa- 
rameters for the second local (or global) energy mini- 
mum state of the zigzag chains are listed in Table V. 
Note that the equilibrium structural parameters listed 
in Table IV were obtained by unrestricted structural re- 
laxations using the conjugate gradient method (see Sec. 
II). Our present calculations corroborate some of these 
previous findings. For example, in the metastable Zr, Ir, 
Pt and Au chains, we find the angle 9 to be 29.6°,28.4°, 
25.8°, 28.2°, d x = 4.25, 4.00, 4.27, 4.60 A, and spin mo- 
ment m s — 0.30, 0.28, 0.00, 0.00/is/atom, respectively. 
In these metastable Zr, Ir, Pt and Au zigzag chains, 
the total energy is, respectively, 1.02, 0.21, 0.41, 0.33 
eV/atom higher than the ground state triangular zigzag 
chains. However, we don't find a second local (or global) 
energy minimum state in the Ru and Os zigzag chains, in 
contrast to the previous studies i 16 i 17 The discrepancy on 
the Ru chain between the present and previous^ 7 - studies 
could be attributed to the fact that highly accurate PAW 
potential rather than norm-conserving pseudopotential, 
is used here, while the difference on the Os chain might 
be caused by the use of the faster but less accurate norm- 
conserving pseudopptential linear combination of atomic 
orbitals method in Ref. 16. Surprisingly, the second en- 
ergy minimum state in the Nb, Mo, Tc, W and Re chains 
(Table V) is in fact the global energy minimum, i.e., its 
total energy is below the corresponding minimum energy 
listed in Table IV. This result for the W chain is in agree- 
ment with the previous study of Ref. Il6l Moreover, the 
bondlength d\ of this second minimum state is not much 
larger than that of the first minimum state and the angle 
a is not close to 30.0°. These results appear to be consis- 
tent with the observation that only Ir, Pt and Au could 
form an atomic chain in the break-junction experiments. 



TABLE VI: Spin (m s ) and orbital (m ) magnetic moments (in 
/is /atom) of the magnetic 4d and 5d transition metal zigzag 
chains in the equilibrium structures (Table IV) with magne- 
tization parallel (m|| z) and perpendicular (m|| x, m|| y) 
(see Fig. 1) to the chain axis from fully relativistic charge 
self-consistent calculations. Superscript a denotes the orbital 
moments on two neighboring atoms are antiparallel, though 
the system is in spin ferromagnetic state. Superscript y means 
that the orbital moment is along the j/-axis, though the spin 
moment is along the a;-axis, i.e., the spin and orbital moments 
are noncollinear. 
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Ru 
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1 
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0.151 
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1 
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(AF) 
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0.263 


0.024 
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0.009 


Rh 


(FM) 


1 


356 


0.338 


1.321 


0.205 


1 
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0.085 


Pd 


(FM) 








0.182 


0.045 
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0.022 




(AF) 





226 


0.018 


0.237 


0.070 





229 


0.049 










5d metals 










Ta 


(AF) 















111 


-0.266 


W 


(FM) 





261 


-0.042 
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0.002 


Re 


(FM) 





517 


-0.018 


0.506 
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0.145 


Os 


(AF) 
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-0.160 


Ir 


(FM) 
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0.433 
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0.458 


Pt 


(AF) 















139 


0.325 


Au 


(AF) 
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0.135 



It could be worthwhile to search for the atomic chains in 
the break-junction experiments using Zr. Finally, the Ir 
and Pt zigzag chains were reported to have a high-spin to 
low-spin transition near the local energy minimum . 16 ! 55 
In the present studies, the spin magnetic moment for Ir 
(Zr) in the ladder-like structure is 0.629 (0.162) /is/atom 
but becomes 0.285 (0.295) /is/atom in the elongated en- 
ergy minimum structure. In the W and Re zigzag chains, 
we found a magnetic to nonmagnetic transition from the 
first energy minimum to the second energy minimum. 
For the Nb, Mo, Tc, Pt and Au zigzag chains, both first 
and second energy minimum states are nonmagnetic, and 
therefore, no high-spin to low-spin transition occurs. 

When the SOC is taken into account, not only the spin 
magnetic moments would depend on the magnetization 
direction but also the orbital magnetic moments would 
appear. In the 3d TM chains, the spin magnetic moments 
are hardly affected by the SOG22 because of the smallness 
of the SOC in these systems. In contrast, in the 5d TM 
chains, the SOC is so large that it not only would affect 
the size of the magnetic moments but also could suppress 
or induce magnetism itself, depending on the magnetiza- 
tion orientation, as mentioned already in Sec. Illb. The 
magnetic moments in the magnetic zigzag chains for the 
magnetization along three coordinate axes from fully rel- 
ativistic charge self-consistent calculations are listed in 
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FIG. 8: Scalar-relativistic (a-c) and fully relativistic (d-f) 
band structures of the Re, W, and Ir zigzag atomic chains 
in the FM state. In (d-f), the spin magnetization is along the 
chain direction (i.e., the 2-axis). The Fermi level (the dotted 
horizontal line) is at the zero energy. 
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FIG. 9: (color online) The cohesive energy of the Ad TM 
zigzag chains in the NM, FM and AF states. For comparison, 
the ground state cohesive energy of the corresponding linear 
chains is also plotted (solid circles). The ground state mag- 
netic configuration of the linear chains is labelled as NM or 
FM or AF near each solid circle. 



Table V. We notice that all 5d TM chains exhibit the 
remarkable CMA4^ behavior. Even two Ad (AF Tc and 
FM Pd) zigzag chains show the CMA too. In particular, 
in the FM Ir chain, the magnetism occurs only when the 
magnetization is along the chain. In contrast, in the AF 
Ir chain, the magnetism appears only when the magne- 
tization is parallel to the y-axis (Fig. lb). The orbital 
magnetic moments in the Zr, Tc, Ru, Rh, Ta, Re, Os, 
Ir, Pt, Au chains are rather significant (> 0.1 /is/atom) 
(Table V). In the Rh, Ir and Pt zigzag chains, the orbital 
magnetic moments can be as large as 0.3 [is /atom. As 
in the linear chains (Table II), the orbital moments in 
the zigzag chains depend strongly on the magnetization 
orientation (Table VI). 



B. Band structures and density of states 

The band structures of the FM Re, W and Ir zigzag 
chains are displayed in Fig. 8, as representatives. Com- 
pared with the corresponding band structures of the lin- 
ear chains (e.g., Ir in Fig. 5), the number of bands be- 
come doubled in the zigzag chains because of the doubling 
of the number of atoms. Furthermore, unlike the linear 
chains where the d xy (d xz ) and d x 2_ y 2[d yz ) bands (Fig. 
5d) are degenerate because of rotational invariance, the 
d xy (d xz ) and d x 2_ y 2(d yz ) bands are now split because of 
the strong anisotropy in the x — y plane perpendicular 
to the chain axis. It is clear that the energy bands are 
also highly spin-split and the separation of the spin-up 
and spin-down bands may be correlated with the spin 
magnetic moment. 

As for the linear chains, we calculate the spin- 



polarization (P) and count the numbers of spin-up and 
spin-down conduction bands at the Fermi level in the FM 
zigzag chains, as listed in Table III. The P in the consid- 
ered zigzag chains generally gets reduced when compared 
with that in the linear chains (Table III). Nevertheless, 
the P of the Rh zigzag chain is as large as 0.53. Inter- 
estingly, the sign of the P in the Zr chain changes from 
positive to negative when it transforms from the linear 
to zigzag chain structure. 



V. STABILITY OF LINEAR CHAIN 
STRUCTURES 

Let us now examine the relative stability of the linear 
and zigzag chains by comparison of the total energies of 
the two structures. The ground state cohesive energy of 
the linear chains and the cohesive energies of the zigzag 
chains in the NM, FM and AF states are displayed in Fig. 
9 (for Ad TM) and Fig. 10. (for 5d TM). The cohesive 
energy (E c ) is defined as the difference in the total en- 
ergy between the free atom (E a ) and the chain (E t ), i.e. 
E c = E a — E t . A positive value of the E c means that the 
formation of the chain from the free atoms would save 
energy, i.e., the chain would be stable against breaking 
up into free atoms. The total energies of the free atoms 
are calculated by the cubic box supercell approach with 
the cell size of 15 A. The electronic configurations used 
for Ad TM are Ad^s 2 (Y), 4d 3 5s 1 (Zr), Ad 4 5s 1 (Nb), 
4d 5 5s 1 (Mo), Ad 6 5s 1 (Tc), 4d 7 5s x (Ru), 4d 8 5s 1 (Rh), 
Ad 9 5s 1 (Pd) and 4d 10 5s 1 (Ag). And for 5d TM are 5d 1 6s 2 
(La), 5d 3 6s 1 (Hf), bd^Qs 1 (Ta), hd^&s 1 (W), bd^Qs 1 (Re), 
SdHs 1 (Os), 5d 8 6s 1 (Ir), 5d 9 6s 1 (Pt) and 5d 10 6s 1 (Au). 
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FIG. 10: (color online) The cohesive energy of the 5d TM 
zigzag chains in the NM, FM and AF states. For comparison, 
the ground state cohesive energy of the corresponding linear 
chains is also plotted (solid circles). The ground state mag- 
netic configuration of the linear chains is labelled as NM or 
FM or AF near each solid circle. 



We note that in all cases, the ground state cohesive en- 
ergy of the linear chain is smaller than that of the zigzag 
chain (Figs. 9 and 10). This suggests that the Ad and 5d 
linear chains are unstable against the zigzag structural 
distortion, as may be expected from the Peierls instabil- 
ity of linear one-dimensional monoatomic metals4£ The 
difference in the ground state energy between the lin- 
ear and zigzag structures for all the 4d and 5d elements 
is rather large, ranging from 0.8 to 2.0 eV/atom. This 
shows that the free standing 4d and 5d TM linear chains 
would not be the stable state, and the linear chains may 
occur only in constrained conditions such as on the steps 
on a vicinal surface^ and under tensile stress in the break- 
point experiment s 8 ' 49 " — . Interestingly, a recent ab initio 
study£2 showed that alloying the gold nanowires with ce- 
sium could make linear monoatomic chains stable. 



VI. MAGNETIC ANISOTROPY ENERGY 

The total energy as a function of the magnetization 
orientation (6, 4>) of a ID wire may be written, in the 
lowest non- vanishing terms, as 



E t =E Q + sin 2 9 {Ei - E 2 cos 2 



(2) 



where 6 is the polar angle of the magnetization away from 
the chain axis (z-axis) and (f> is the azimuthal angle in the 
x — y plane perpendicular to the wire, measured from the 
x axis. For the free standing linear atomic chains, the 
azimuthal anisotropy energy constant E 2 is zero. The 
axial anisotropy energy constant E\ is then given by the 
total energy difference between the magnetization along 
the y(x) and z axes, i.e., E x = EV - E z (E x = E»). A 



TABLE VII: Total (El), electronic (Ef) and dipolar (Ef) 
magnetic anisotropy energies (in meV/atom) of the 4d and 
5d transition metal linear chains. If E\ is positive, the easy 
magnetization axis is along the chain; otherwise, the easy 
magnetization axis is perpendicular to the chain. 
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0.141 


Tc 








7.228 


7.186 
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Ru 


-11.99 
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Rh 
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Pd 
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0.012 














5d metals 








Hf 


0.825 


0.825 


0.000 








W 








-5.235 


-5.283 


0.048 


Re 








-59.94 


-60.01 


0.070 


Ir 


-11.13 
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positive value of E\ means that the chain (z) axis is the 
easy magnetization axis. For the zigzag chains which are 
in the x — z plane, E 2 is not zero and can be calculated 
as the total energy difference between the magnetization 
along the x and y axes, i.e., E 2 = E y — E x . 

The magnetic anisotropy energy for a magnetic solid 
consists of two contributions. One comes from the mag- 
netocrystalline anisotropy in the electronic band struc- 
ture caused by the simultaneous occurrence of the elec- 
tron spin-orbit interaction and spin-polarization in the 
magnetic system, and ab initio calculation of this part has 
already been described in Sec. II. The other is the magne- 
tostatic (or shape) anisotropy energy due to the magnetic 
dipolar interaction in the solid. The shape anisotropy 
energy is zero for the cubic systems such as bcc Fe and 
fee Ni, and also negligibly small for weakly anisotropic 
solids such as hep Co. However, for the highly anisotropic 
structures such as magnetic Fe and Co monolayers ) 5 ^ 54 
the shape anisotropy energy can be comparable to the 
electronic MAE, and therefore cannot be neglected. For 
the collinear magnetic systems (i.e. m q / /m q > ), this mag- 
netic dipolar energy Ed is given by (in atomic Rydberg 
units)£2 



Ed = E 



(3) 



where M qq > is called the magnetic dipolar Madelung con- 
stant which is evaluated by Ewald's lattice summation 
technique^. The speed of light c = 274.072, and m q is 
the atomic magnetic moment on site q in the unit cell. 
Note that in atomic Rydberg units, one Bohr magneton 
(/is) is y/2. Therefore, as noted recently in Ref. HH, 
the Ed for the multilayers obtained previously by Guo et 
aZj 53 i 54 is too small by a factor of 2. 

The calculated E d 's for the linear and zigzag chains are 
listed in Tables VII and VIII, respectively. Tables VII 
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and VIII show that in both the linear and zigzag chains 
and in both the FM and AF states, the E d, s are much 
smaller than the electronic contributions (E e ), being in 
strong contrast to the case of the 3d TM chains^ This 
is because the magnetization here is significantly lower 
and the equilibrium bond length becomes larger, com- 
pared with that of the 3d TM chains^ Furthermore, 
they always prefer the chain direction (z axis) as the 
easy magnetization axis. Therefore, any perpendicular 
magnetic anisotropy must originate from the electronic 
magnetocrystalline anisotropy. 

The calculated E e, s of the linear and zigzag atomic 
chains are also listed in Tables VII and VIII, respectively. 
Table VII shows that in the FM linear chains at equilib- 
rium, the E e would favor a perpendicular anisotropy in 
the Zr, Ru, Pd and Ir chains but prefer the chain axis 
in the Rh and Hf chains. In the AF state, in contrast, 
the Mo, W and Re linear chains would have the easy 
axis perpendicular to the chain while only the Tc linear 
chain prefer the axial anisotropy. Remarkably, the FM 
Ru, Rh and Ir as well as AF Tc and Re linear chains 
have a large total anisotropy energy (£?*) (see Table VII) 
of ~10 meV/atom. In particular, the E l of the AF Re 
linear chain is as large as -60 meV/atom. Ab initio cal- 
culations of the E e of the 4c? TM linear chains have been 
reported recently^, and our present results for the equi- 
librium bondlengths (Table VII) agree rather well with 
these earlier calculations (Fig. 1 in Ref. l44h . 

The electronic anisotropy energy for the selected linear 
Ad and 5d atomic chains is displayed as a function of 
bond length in Fig. 11. It is clear that in several selected 
linear chains, the magnitude of the E e generally increases 
with the bondlength (Fig. 11a), like the spin and orbital 
magnetic moments (Fig. 4). For example, the E e of the 
Rh chain increases from 7.0 me V/ atom at the equilibrium 
bondlength (2.25 A) to 37.3 meV/atom at bondlength 
of 3.0 A. Several chains also undergo interesting spin- 
reorientation transition as the bondlength is elongated. 

When elongated, for example, the FM Zr, Nb, Ru, 
and Ir linear chain would undergo a spin reorientation 
transition from the perpendicular to along the axial di- 
rection at the bondlength of ~2.75A, -2.82 A, ~2.65A, 
and ^2.45A, respectively. In contrast, the magnetiza- 
tion of the Ta chain transits from that along the axis 
to the perpendicular direction at ^2.85A. Furthermore, 
many elongated chains have a gigantic anisotropy energy 
of ~20 meV/atom (Fig. 11). 

Table VIII shows that the size of the axial anisotropy 
energy {E\) in the zigzag chains is large and is generally 
comparable to that in the linear chains (Table VII) . How- 
ever, unlike the linear chains, there is also the pronounced 
anisotropy (E2) in the x — y plane perpendicular to the 
chain axis in many zigzag structures (Table VIII) . In the 
FM Y, AF Ru and Os zigzag chains, the magnitude of the 
E l 2 is even larger than that of E\. In the FM Y, Mo, AF 
Ru, FM Rh, FM Re, and FM Ir zigzag chains, the easy 
axis is along the chain direction. In the FM Ru as well as 
FM and AF Pd chains, the easy axis is perpendicular to 



the zigzag plane. In the Zr, W, and Os chains, the easy 
axis is in the zigzag plane but perpendicular to the chain 
axis. Ab initio calculations for only the Ir and Pt zigzag 
chains have recently been reported*^ However, in Ref. 
[iH the easy axis is reported to be along the x-axis. This 
discrepancy may be due to the pronounced difference in 
the equilibrium zigzag structure between the present and 
previous calculations. Furthermore, we find the Pt zigzag 
chain to nonmagnetic. 

Structural transformation from the linear to zigzag 
structure has profound effect on magnetism in the Ad 
and 5d TM nanowires. This transformation not only in- 
duces (or suppresses) magnetization in, e.g., the Y and 
Os chains (the Tc and Hf chains), as mentioned already 
in Sec. IV, but also causes spin reorientation transition 
in, e.g., the Re and Ir chains (Tables VII and VIII). Note 
that the linear AF Re chain has a gigantic perpendic- 
ular anisotropy energy of -60.0 meV/atom (Table VII). 
However, upon transition to the zigzag structure, the AF 
state disappears, and, instead, the FM state appears with 
the magnetization switched to be along the chain axis. 



VII. CONCLUSIONS 

We have performed an extensive ab initio study of 
the physical properties of both linear and zigzag atomic 
chains of all Ad and 5d transition metals within the GGA 
by using the accurate PAW method. First, the atomic 
structures were determined. All the TM linear chains are 
found to be unstable against the corresponding zigzag 
structures. All the TM chains except Nb, Ag and La, 
have a stable (or metastable) magnetic state in either 
the linear or zigzag or both structures. Magnetic states 
appear also in the Nb and La linear chains when the 
chains are sufficiently elongated. The spin magnetic mo- 
ments in the Mo, Tc, Ru, Rh, W, Re chains could be 
large (>1.0 /its /atom). Structural transformation from 
the linear to zigzag chains can suppress the magnetism 
already in the linear chain, induce the magnetism in the 
zigzag structure, and also cause a change of the magnetic 
state (FM to AF or vice verse). 

With the SOC included, our calculations show that 
the orbital moments in the Zr, Tc, Ru, Rh, Pd, Hf, 
Ta, W, Re, Os, Ir and Pt chains could be rather large 
(>0.1 /^s/atom). Importantly, large magnetic anisotropy 
energy (>1.0 meV/atom) is found in most of the mag- 
netic TM chains, suggesting that these nanowires could 
have important applications in ultrahigh density mag- 
netic memories and hard disks. In particular, giant mag- 
netic anisotropy energy (>10.0 meV/atom) could appear 
in the Ru, Re, Rh, and Ir chains. Furthermore, the mag- 
netic anisotropy energy in several linear chains could be 
as large as 40.0 meV/atom when the chains are under 
sufficiently large tensile strain. A spin-reorientation tran- 
sition occurs in the Ru, Ir, Ta, Zr, La and Zr, Ru, La, 
Ta and Ir linear chains when they are elongated. Re- 
markably, all the 5c? as well as Tc and Pd chains show 
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TABLE VIII: The total (E{, E$), electronic {El , _E|) and dipolar (E d , E%) magnetic anisotropy energy constants (in meV/atom) 
as well as the easy magnetization axis (M) of the 4d and 5d transition metal zigzag chains. Ei — E v - E z ; E2 — E v - E x , see 
Equ. (2). 
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FIG. 11: (color online) Magnetocrystalline anisotropy energy 
(El) of the selected 4d and 5d transition metal linear atomic 
chain as a function of interatomic distance. The upper panels 
contains the TM linear chain with larger MAE. A positive 
value of El means that the magnetization would be parallel 
to the chain axis whilst a negative value would means that the 
easy magnetization axis would be perpendicular to the chain. 



the fascinating behavior of the so-called colossal magnetic 
anisotropy. 45 Finally, the electronic band structure and 
density of states of the nanowires have also been calcu- 
lated mainly in order to understand the electronic origin 
of the large magnetic anisotropy and orbital magnetic 
moment as well as to calculate the conduction electron 
spin polarization. 
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